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Current computational capability allows a large-scale first-principles calculation, however, the 
heavily folded electronic band structure due to a big supercell suppresses an easy visualization of the 
effect of symmetry breakers introduced to the reference system. While an unfolding method has been 
proposed, a general formula is not yet available for a basis set of non-orthogonal atomic-like orbitals 
that is ideal for the unfolding concept due to the translational symmetric nature. In this study, we 
derive a simple formula applicable for both orthogonal and non-orthogonal localized orbitals. Codes 
using a basis set of linear combination of atomic-like orbitals, where the order-JV method can be 
implemented much easier, can benefit from this approach in unfolding band structures. This method 
is applied for analysis of surface state with ZrB2(0001) surface and provides a way to examine how 
good a surface-sensitive measurement, like ARPES, can represent the bulk property. The effect of 
vacancies at various concentration levels in ZrB2 will also be discussed. 

PACS numbers: 71.15.-m,71.20.-b,79.60.-i 



I. INTRODUCTION 

The success of Kohn-Sham framework of density func- 
tional theory (DFTjB has led to extensive realization 
of solid state physics from the electronic band struc- 
tures. One of the most utilized strategy to perform a 
first-principles calculation is the adoption of a supercell 
approach that makes studies of various forms of imper- 
fections possible, for example, impurities, vacancies, a 
charge density wave, a spin density wave, and/or an or- 
bital ordering introduced to extended systemsP Not re- 
stricted by the periodic boundary condition, the Bloch 
theorem can also be applied to study non-periodic sys- 
tems, thank to the feasible computation of a big enough 
supercell that also allows one to investigate a system 
highly breaking the translational symmetry such as a 
presence of a surface. 3 However, a bigger supercell im- 
plies a smaller 1st Brillouin zone (BZ), as shown in Fig.JT] 
that highly suppresses an easy visualization of the effect 
of symmetry breakers introduced to the reference system 
due to the heavily folded bands in the 1st BZ that no 
longer resemble the original band structure anymore. A 
pioneering concept has been proposed to represent the 
spectral function, the imaginary part of retarded one- 
particle Green function, in the basis of the reference sys- 
tem, not the one of the supercell, allowing a direct com- 
parison between first-principles calculations and angle re- 
solved photoemission spectroscopy (ARPES) P 

While the unfolding method has been introduced via 
the use of informative Wannier functions and gives rise 
to many detailed insights of physical understanding on 
various systems)^' the construction of Wannier func- 
tions needs an additional effort from resulting Bloch 
states and might cost a non-negligible time for large- 
scale systems. Besides, the "unfolding" would physically 




FIG. 1: Real-space unit cells of (a) normal-cell ZrB2, (b) 
(4x4x4) supercell with a B vacancy, and (c) 13-layer slab. 
The corresponding 1st Brillouin zones are presented in (d), 
(e), and (f), respectively. Note that a nearly the same amount 
of states are accommodated in all illustrated zones, indicating 
a heavily folded feature of the supercell band structures. 



make sense if the basis functions are almost identical be- 
tween different reference unit cells, therefore, the Wan- 
nier function should be used with caution in unfolding 
the band structure, owing to the gauge freedom dur- 
ing the construction.^ The unfolding method has also 
been introduced to the plane wave basis set, allowing 
an alternative option with different emphasis on physi- 
cal properties.^ In fact, both implementations imply an 
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orthonormal relationship for the smaller reference sys- 
tem. Given the translational symmetry has been broken 
in the actual supercell due to the presence of the sym- 
metry breaker, it would be interesting to explore the un- 
folding method by using a translational symmetric basis, 
where the orthogonality nature must be brokenP^Hm 

For studying a large-scale system, it is desirable to 
allow a linear scaling (order-iV) computation that can 
highly optimize understanding of large-scale systems. It 
has been developed by using the basis of linear com- 
bination of atomic-like orbitals (LCAO)PHH Although 
the orthonormal relationship introduced for the unfold- 
ing method using Wannier functions cannot be directly 
applied for the LCAO basis without further orthonormal- 
izing the basis functions, it provides another opportunity 
to explicitly take the overlap matrix of atomic-like or- 
bitals, which are ideal for the unfolding concept, into 
account. In order to complement the current computa- 
tional capability, it is timely to introduce a more general 
formula to optimize a deeper understanding of the effect 
of symmetry breakers in the band structures of large- 
scale systems with the LCAO basis. 



II. OVERVIEW AND COMPUTATIONAL 
DETAIL 

In this study, we derive a general formula in detail for 
unfolding electronic band structures applicable for both 
orthogonal and non-orthogonal basis sets via an intro- 
duction of the overlap matrix of localized orbitals. From 
that, it becomes clear on how an orthonormal basis set 
can simplify the calculation and on how good an approxi- 
mation is if one chooses to neglect the non-orthogonality 
nature with a non-orthogonal basis set in order to re- 
duce the computational time, which is illustrated with 
a study of (2x2x4) supercell without introducing any 
symmetry breaking. This method is applied for study- 
ing the proper spectral weight of ZrB2 with the presence 
of (lxl) and (2x2) unit cells of (0001) surface. The 
revealed spectral weight allows us to study if a surface 
component can represent the bulk band structure. Ad- 
ditionally, we have also investigated the effect of a low- 
concentration vacancy by introducing one Zr or one B 
vacancy in Zr 64 B 12 8 (Zro.9s437.5B2 and ZrBj. 954375) and 
have examined the applicability of a rigid band shift ap- 
proximation upon doping for a metallic system. 

The DFT calculations within a generalized gradi- 
ent appro ximation^! are performed by the OpenMX 
codej 12 ! 21 ! which is based on norm-conserving pseudopo- 
tentials generated with multi reference energies and opti- 
mized pseudoatomic basis functions. The cut-off energy 
is set to 270 Ry and the basis functions of Zr— s3p3d2 and 
B— s4p2dl are chosen with cut-off radius of 7 Bohr. The 
lattice constants of bulk ZrB 2 converge to a=3.174A and 
c=3.55A with a 8x8x5 mesh of k points. The structure 
can be found in Fig. [T] The k points of supercell cal- 
culation are chosen to have the same quality of k mesh 



as the bulk one. The atomic positions are relaxed and 
the residual force on each atom is less than or equal to 
3 x 10~ 4 Hartrcc/Bohr. The bulk lattice constants are 
in a good agreement with experimental dataP^l 



III. DERIVATION OF PROPER SPECTRAL 
WEIGHT IN A NON-ORTHOGONAL BASIS 

We now introduce a simple formula to unfold the band 
structure with the LCAO basis. The idea is to represent 
the spectral function in the LCAO basis of the reference 
system where a smaller normal unit cell is adopted. It is 
obvious that the unfolding process would reproduce the 
same result of the calculation of normal cell if the sym- 
metry breaking is switched off in the supercell. There- 
fore, the strength of each band's coupling to the symme- 
try breaker can be observed via the calculated spectral 
weight. One essential property of the LCAO basis is the 
non-orthogonality nature whose strength can be found 
from the overlap matrix elements of localized orbitalsP^ 
Since the overlap matrix elements and the LCAO coeffi- 
cients of eigenstates at desired K points can be obtained 
from a regular band structure calculation of a big su- 
percell, like the one provided by the OpenMX code, we 
assume those are available information and will derive a 
formula directly containing those two terms. 

One can choose the K path to be outside the 1st BZ of 
the supercell to allow a comparison with the band struc- 
ture along the k path in the normal-cell BZ. Although 
the K point outside the 1st BZ gives the same eigenstates 
(therefore, the same LCAO coefficients) as the one inside 
the 1st BZ connected by reciprocal lattice vectors (G) of 
the supercell (will be referred to 8}. ,[k])> the correspond- 
ing spectral weight is modified due to the new represen- 
tation. The proper spectral weight, A(lo), represented 
by the conceptual normal-cell eigenstate \kj) (lowercase 
letters) can be obtained via this formula: 

J2(kj\A(u)\kj) 
kj 

= Y, (kj\km)S^ n (k)(kn\A\km')S m } nl (k)(kn'\kj) 

kjmnm' n' 

= J2 (kn^km)S- 1 n (k)(kn\A\km')S-} n ,(k) 

kmnm'n' 

= E S m ) n ,(k)S n , m {k)S^ n {k){kn\A\km') 

kmnm'n' 

= J2 S ™l(k)(kn\A\km). 

kmn 

The overlap matrix, S(k), does not follow the or- 
thonormal relationship in general and can not be ne- 
glected. By inserting known information, such as 
a supercell representation (uppercase letters) in real 
space (Y,rmr'M'\ RM ) S rm,r>M>( R ' M '\) or the super- 
cell eigenstate Q2kj \KJ){KJ\) in this study, the spec- 
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tral function at each k becomes: 
j 

= E S m i{k){kn\KJ){KJ\A\KJ){KJ\km) 

mnK J 

= E S-i(k){kn\KN}C§ J (KM\kTn}C§ J *A KJ<K ^ 

mnK J M N 

(2) 

C§ J denotes the LCAO coefficients of an eigenstate 
\KJ), such that \KJ) = J2 N C§ J \KN). 

Following the concept of unfolding, it is desirable to 
relabel the supercell orbital \RN) by the normal-cell \rn) 
and find the corresponding inner product: 

(rn\RN) 

= (m\R + r (N),n'(N)) 

= E(rn|fc'n)(fc'n|fcV(iV))(fc'n'(^)|i? + ro(A r ),"'(A r )) 



localized orbitals. 

The spectral weight of an eigenstate of the supercell, 
Akj,kj(u), is just a delta function at the correspond- 
ing eigenvalue, therefore, one can calculate the spectral 
weight individually, for example, only the eigenstates 
within an energy window is needed in actual calculation. 
L(l) denotes the number of K(k) points in the supercell 
^(normal-cell) first Brillouin zone and j is the ratio of the 
volume of 1st BZ of the supercell to the normal-cell one. 
In the last step of derivation in Eqn. [5] the overlap ma- 
trix element (R' M\r'n' (N)) has been reduced to the one 
with R' = for a periodic system: 

E e- lKR 'e tkr '(R'M\r'n'(N)) 



E 



o ik'(r-R-r (N)) 



I 



-S n n'{N)(k'). 



(3) 

N labels the supercell orbital and n(N) indicates which 
orbital the N belongs to in the normal cell. In order to 
employ Eqn. [3j we can insert real-space lattice vectors 
in Eqn. [2j Finally, the spectral function in the normal- 
cell representation can be derived in Eqn. [5] by directly 
containing LCAO coefficients and the overlap matrix of 



r'R' 



-i(k-G)R'ikr' I pi MW'r,' 



(R'M\r'n'(N)) 



r'R' 

E e iGR ' E e- lk ^'-^ (R'M\r'n'(N)) (4) 

R> r' 

E e tGR ' E eifcr ' (OM|r'n'(A0) 

R' r' 

L J2e lkr '(0M\r'n'(N)). 



S(R = 0, M; r, n(N)) denotes the overlap matrix ele- 
ments including all non-zero contribution after a rela- 
beling with the lattice vector r and the orbital n of the 
normal cell. 



E^fcjM = E S~l(k)(kn\rn)(rn\RN)(RN\KN)C§ J (KM\R'M)(R'M\r'm)(r'm\km^ 



mnKJMNrr'RR' 
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5 HK] C^ J {R'M\r'n'{N))C^ J A KJ , KJ (u) 
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KJMNr' 
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KJMNr 



Clearly, the effect of translational symmetry breaking has been not only recorded in the supercell eigenstates 
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but also connected via the overlap matrix elements that 
are essential to deliver a complete spectral weight in the 
new representation. Note that this formula is in fact 
a more general formula for a basis set of localized or- 
bitals and can be applied for both orthogonal and non- 
orthogonal orbitals by providing the corresponding S. It 
is obvious that the summation can be reduced by an 
adoption of orthogonal basis. Instead of carrying the 
relative phase from the distance between a lattice vec- 
tor r and a r'(N) weighted by the overlapping between 
\0M) and \rn(N)) orbitals to give rise to the correspond- 
ing contribution of \KM) and \KN), only the absolute r 
vector in a supercell is needed. 

Since S decays rapidly due to the localized charac- 
ter of LCAO basis, it could be a good approximation 
to assume an orthonormal relationship, especially for 
occupied bands. However, it should be checked care- 
fully in case of existence of strong effects from the non- 
orthogonality nature in studied cases. Nevertheless, the 
complete proper spectral weight in non-orthogonal ba- 
sis only involves a simple summation with appropriate 
phases and directly uses the information from a regular 
band structure calculation, allowing an easy implementa- 
tion. Finally, an arrangement of this formula by separat- 
ing the summation of orbital index, e.g. • Akj,kj(w) — 
Wm, further provides the corresponding orbital con- 
tribution Wm in the unfolded band structure. 



IV. ASSUMPTION OF ORTHONORMAL 
RELATIONSHIP IN LCAO BASIS 

One can study how good an approximation is by as- 
suming an orthonormal relationship even with the use of 
non-orthogonal LCAO basis. We illustrate an example 
of a (2x2x4) supercell calculation without introducing 
any symmetry breaking. Therefore, the unfolding pro- 
cess is expected to unfold the heavily folded band struc- 
ture back to the normal-cell one exactly. We compare the 
unfolded band structures with and without the assump- 
tion of orthonormal relationship by tuning the overlap 
matrix S in Fig(2] As can be seen clearly, the unoccu- 
pied bands, where the eigenstates deviate from atomic 
orbitals much more, need large values of coefficients on 
non-orthogonal atomic orbitals and show strong disagree- 
ment. Without symmetry breaking, the deviation of oc- 
cupied bands is small and the shadow bands show negli- 
gible intensity. Therefore, it could be a good approxima- 
tion to assume the orthonormal relationship in unfold- 
ing occupied bands. This approximation is useful for a 
linear scaling calculation, where the one-particle Green 
function is usually obtained, not the LCAO coefficients, 
in estimating the unfolded band structure. Obviously, in 
order to deliver a complete spectral weight, one should 
not neglect the overlap matrix and one should check the 
approximation with caution, especially for the case with 
strong symmetry breaking. 




r m k r 



FIG. 2: Band structure of (2x2x4) supercell calculation of 
ZrB2 without symmetry breaking. The result is presented 
in the first Brillouin zone of (lxlxl) normal cell. Unfolded 
band structures with and without the assumption of orthonor- 
mal relationship are presented by red and black circles, re- 
spectively, and the radius reflects the intensity of the proper 
spectral weight of each supercell eigenstate. The radius of red 
circle reflects how good the approximation is by comparing to 
the reference black circle. The degenerate eigenstates (energy 
difference less than 0.001 eV) are presented by one circle. The 
green dots show the band structure without unfolding from 
the (2x2x4) supercell calculation for comparison. 



V. SURFACE STATE WITH ZIRCONIUM 
DIBORIDE (0001) SURFACE 

A long standing issue on how good a surface-sensitive 
measurement, such as ARPES, can represent the bulk 
properties can be explicitly studied via unfolding. With 
the computational capability nowadays, the surface and 
bulk properties can be revealed via a slab calculation. 
Once the slab is thick enough, the middle part of the 
slab behaves as a bulk even with the presence of the sur- 
face. Therefore, the electronic band structure of such a 
slab exhibits both surface and bulk contributions simul- 
taneously. Of course, a thicker slab is needed for a sys- 
tem with a strong interaction between surface and bulk 
layers. On the other hand, a bigger unit cell of studied 
surface is usually adopted in a common practice due to a 
spontaneous recon struction^ or an interaction with new 
added layers a simple projection of spectral weight 
on terminated layers can only reflect the spectral weight 
on the surface but cannot resolve the folding effect. The 
concept we have introduced is exactly aimed to unfold 
the band structure in relieving the complexity. 

A graphene counterpart, silicene, has recently been 
epitaxially grown on the (2x2) unit cell of ZrB2(0001) 
surface. 26 Without silicene, no surface reconstruction is 
found and only (l xl) u nit cell is needed for studying 
ZrB 2 (0001) surface! 2 ^ Therefore, we use it as a case 



FIG. 3: Electronic band structures of Zr-terminated slabs 
for different layers with the presence of (2x2) unit cell of 
ZrB2(0001) surface. The proper spectral weight of each su- 
percell eigenstate in bulk representation is proportional to the 
radius of the corresponding circle. For a better visualization, 
the degenerate eigenstates (energy difference less than 0.001 
eV) are presented by one circle. The dispersion curve indi- 
cated by the red arrow can be identified as the surface state 
and the gap opening deviating from the bulk band structure 
is pointed out by the green arrow. 

study to investigate the proper spectral weight by per- 
forming a slab calculation with the (2x2) unit cell of 
ZrB 2 (0001) surface and unfold the band structure to 
the (lxlxl) bulk representation. One can expect the 
unfolding method can recover the result of (lxl) unit 
cell of the surface from the introduced redundant (2x2) 
unit ce ll and can have an easy comparison with available 
results! 2112 ^ 

After relaxing atomic positions in a slab with an addi- 
tional empty layer (composed of LCAO basis of B atom 
with a flat zero potential) added into Zr terminated sur- 
face to have more complete information in the vacuum, 
band structures with different number of layers are pre- 
sented in Fig. [3j Clearly, a band that cannot been iden- 
tified in the normal-cell calculation [c.f. Fig. [5] crosses 
the Fermi energy between T and M in four different slab 
calculations and can be attributed as a surface state dis- 
persion (red arrow). Note that the comparison of the 
band structure to the normal-cell one becomes easier in 
identifying the surface state due to unfolding the 1st BZ 
of (2x2) surface unit cell to the bulk one. Several shadow 
bands and gap opening can be observed due to the pres- 
ence of the surface. Interestingly, while the surface state 
is converged with the number of layers and is good for a 
surface study, a gap opening between M and K at the 
Fermi energy (green arrow) is present even with the 19- 
layer slab calculation that deviates from the bulk band 
structure. One can see the gap closing converges slowly 
in Fig. [3j indicating a requirement of much more layers to 



FIG. 4: Electronic band structure of Zr-terminated 59-layer 
slab with the presence of (lxl) unit cell of ZrB2(0001) sur- 
face. The radii of circles are proportional to the positive con- 
tribution of proper spectral weight in the bulk representation. 
The red circles indicate the contribution of the two outermost 
Zr layers (four for two surfaces in a slab). The blue circles 
show the weight from the two outermost B layers (four for two 
surfaces in a slab). The rest contribution is presented by the 
green circle and reflects the bulk property. For a better vi- 
sualization, the degenerate eigenstates (energy difference less 
than 0.001 eV) are presented by one circle and the radii of 
red and blue circles are doubled for comparison. 



capture the bulk band structure, highlighting the needs 
of a larger scale calculation and the unfolding method for 
it. It also reflects a deeper probe by a measurement is 
needed to observe a crossing dispersion between M and 
K below the Fermi energy. 

An advantage of the LCAO basis is the real-space rep- 
resentation that allows a useful analysis on the layer con- 
tribution of a slab in the band structure. In order to 
estimate the spectral weight of the surface, the positive 
contribution of Wm from the outermost layers is studied 
for a (lxl) unit cell of (0001) surface with 59 layers (30 
Zr layers and 29 B layers). The one from two outermost 
Zr layers is shown in Fig. [4] by red circles. The outer- 
most Zr layers not only contribute the spectral weight on 
the band of the surface state but also partly follow the 
normal-cell band structure. With unfolding, it is easy 
to see where the red circles deviate from the bulk bands 
(green circles) around the Fermi energy. Blue circles show 
the contribution from two outermost B layers. It is clear 
that the outermost ZrB 2 layers have already participated 
the bulk band structure by falling a weight on the con- 
tribution of rest layers (green circles). This advocates 
even a surface-sensitive measurement, like ARPES ex- 
periment, can access bulk band structures from several 
outermost layers. Of course, how deep the layer can re- 
produce the band structure of the middle part of a slab 
is case dependent and can be checked with a slab calcu- 
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FIG. 5: Electronic band structures of a Zr vacancy intro- 
duced to the indicated supercells of ZrB2. The (lxlxl) 
normal-cell band structure containing one Zr atom and two B 
atoms is presented as a reference for comparison. The proper 
spectral weight of each supercell eigenstate is proportional to 
the radius of the corresponding circle. For a better visualiza- 
tion, the degenerate eigenstates (energy difference less than 
0.001 eV) are presented by one circle. For the (4x4x4) su- 
percell case, although an overall spectral weight converges to 
the normal-cell one, gap openings still remain as indicated by 
the green arrow. 



FIG. 6: Electronic band structures of a B vacancy introduced 
to the indicated supercells of ZrB2. The (lxlxl) normal-cell 
band structure containing one Zr atom and two B atoms is 
presented as a reference for comparison. The proper spectral 
weight of each supercell eigenstate is proportional to the ra- 
dius of the corresponding circle. For a better visualization, the 
degenerate eigenstates (energy difference less than 0.001 eV) 
are presented by one circle. For the (4x4x4) supercell case, 
although an overall spectral weight converges to the normal- 
cell one, a gap opening still remains as indicated by the green 
arrow. 



lation in detail via unfolding. 



VI. VACANCY STATE IN ZIRCONIUM 
DIBORIDE 

A good approximation for studying the doping effect 
on the electronic band structure is the rigid band shift 
one that simply readjusts the Fermi energy of the normal- 
cell calculation according to the amount of electron or 
hole doping. This approximation makes sense if the in- 
troduced dopant interacts weakly with rest of the system. 
The applicability of this approximation can be visual- 
ized by comparing the calculations between a normal-cell 
band structure and a big supercell one via unfolding, not 
just a common comparison of density of states. One can 
expect as the supercell size becomes larger and lager, the 
normal-cell band structure can be recovered with a low- 
concentration disordered symmetry breaker, as can be 



observed by ARPES experiment, not the heavily folded 
bands fully occupying the supercell BZ. In Fig. [5] and |6j 
we show the band structures of ZrB2 with one Zr and one 
B vacancy (replaced by an empty atom) [c.f. Fig. [I], re- 
spectively, in the indicated bigger supercells and compare 
to the (lxlxl) Z1-B2 whose 1st BZ is adopted to repre- 
sent all results. As can be seen, the band structure be- 
comes simpler and simpler with increased supercell sizes. 

Not immune to the ordered vacancy we introduced to 
ZrB2, especially the one in the smaller (2x2x4) super- 
cell, prominent intensity of shadow bands and gap open- 
ings with considerable sizes can be observed, reflecting 
the strength of each band's coupling to the translational 
symmetry breaker, which is the vacancy in this case. The 
band structures of (3x3x3) and (3x3x4) supercells show 
a tendency of convergence toward the normal-cell one. 
However, a flat band around the Fermi energy in Fig. [6] 
can be found that can significantly modify the conduc- 
tivity for insulators due to the originally gapped band 
structure. For a metallic ZrB2, the effect of vacancy 
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states is not profound. Indeed, while the supercell size 
becomes larger, the normal-cell band structure can be re- 
covered without prominent shadow bands and gap open- 
ings, advocating a rigid band shift approximation could 
be a good approximation if an overall band structure is 
concerned in the study of ZrB2 with low-concentration 
vacancies. Nevertheless, one should realize the spectral 
weight does deviate from the normal-cell one more or less 
due to the vacancy even in the lowest concentration level 
we have studied, namely ZrBi. 984375, and could perturb 
several physical properties, for example, a gap opening 
remains as indicated by the green arrow in Fig. [6] due to 
the ordered vacancy. Note that such gap opening can be 
observed easily from the unfolded band structure. 

VII. SUMMARY 

In summary, we offer a general formula to unfold first- 
principles band structures applicable for both orthogonal 
and non-orthogonal localized orbitals as basis functions. 
Codes using a LCAO basis set, which is suitable for a lin- 
ear scaling calculation and for the unfolding concept, can 
directly implement such formula with a simple summa- 
tion from a regular band structure calculation. Different 
from an orthogonal basis, the overlap matrix of local- 



ized orbitals plays an important role in connecting the 
recorded information of translational symmetry break- 
ing in eigenstates. The LCAO basis contains real-space 
information that allows one can easily study the effect 
of symmetry breakers, such as the surface and the va- 
cancy, in real space. By using ZrB2 as a case study, 
we demonstrated unfolded band structures in the repre- 
sentation of the normal-cell 1st BZ with the presence of 
different concentration of vacancies and the surface. A 
slab calculation allows a straightforward way to compare 
with the band structure from a surface-sensitive measure- 
ment. By unfolding the band structure, the comparison 
with bulk states becomes easier in allowing identifying 
the surface state and the contribution of different layers 
can be further studied. In fact, several outermost lay- 
ers are found to be able to exhibit most bulk band dis- 
persion quite well. It is interesting to see if experiment 
can observe a crossing dispersion between M and K be- 
low the Fermi energy that belongs to the bulk property. 
Our method also shows a rigid-band shift approximation 
physically makes sense for a metallic ZrB2 upon doping 
at low-concentration level. Benefiting from the unfolding 
method, a detailed study in the electronic band struc- 
ture of a large-scale system is promising and the unfolded 
band structure can highly complement the information 
one can learn from density of states. 
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